Inversion techniques for optical conductivity data 
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Optical data is encoded with information on tlie microscopic interaction between charge carri- 
ers. For an electron-phonon system, the Eliashberg equations apply and a Kubo formula can be 
used to get the infrared conductivity. The task of extracting the electron-phonon spectral density 
a^F(ijj) from data is rather complicated and, thus, simplified but approximate expressions for the 
conductivity have often been used. We test the accuracy of such simplifications and also discuss the 
advantages and disadvantages of various numerical methods needed in the inversion process. Normal 
and superconducting state are considered as well as boson exchange mechanisms which might be 
appficable to the High- Tc oxides. 
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I. INTRODUCTION 

The interaction Hamiltonian between electrons and 
phonons involves a complicated matrix element or cou- 
pling function gk.k'.u which describes the scattering of 
an electron initially in the state |k) to any final state 
|k') through the exchange of a phonon uj^Ck' — k). Here 

is a phonon branch index and the momentum trans- 
fer k' — k can fall outside the first Brillouin zone and 
so phonon Umklapp processes enter. In a real metal 
the Bloch states of the band structure can be compli- 
cated and this is reflected in the electronic state labeled 
by |k). Fortunately, many important properties of an 
electron-phonon system require for their understanding 
only a Fermi surface to Fermi surface average of the cou- 
pling, namely the functioni 

a^FiLu) = ^ g i?,(k' - k)S (£k - A*) S (£k' - m) , 

where /i is the chemical potential, Ck the electron energy, 
7V(/x) the electron density of states and -Biy(k' — k) the 
phonon spectral function. For example, in the Eliashberg 
formulation'^ of superconductivity based on Migdal's the- 
orem for electron-phonon vertex corrections, it is a^F{uj) 
that enters. For the infrared conductivity another, some- 
what different weighting of 5k,k',iy comes in and the re- 
sulting function of to is usually called the transport spec- 
tral density denoted a'^j.F{u})i^ Here we will not deal di- 
rectly with these differences. An important goal of exper- 
iments in conventional superconductors has been to de- 
termine the electron-phonon spectral function a'^F{uj).'^'^ 
This has been successfully accomplished for a large num- 
ber of conventional materials using tunneling data and 
the inversion technique of McMillan and Rowell.'^ In a 
few cases the infrared optical conductivitji^ was also 
used and excellent agreement with tunneling results was 
found. 

Extensions to the consideration of the A15 compounds 
revealed that additional features of the band structure 



such as the energy dependence of the electronic den- 
sity of states N{e) can also be important. ^'^ More re- 
cently the optical data in the alkali doped Cgo com- 
pounds has been invertediS and found consistent with 
its superconductivity. When experimentally determined 
electron-phonon spectral functions are compared with 
first principle band structure calculations extended to 
include electron-phonon interaction good agreement is 
obtainediSii 

In dealing with the high-Tc oxides several complica- 
tions immediately arise. First, their superconductivity is 
not generally believed to be due to the electron-phonon 
interaction. A consensus exists that the gap has d-wave 
rather than s-wave symmetry and comes out as a re- 
sult of strong correlation effects. A natural explana- 
tion for the d-wave gap is found in the antiferromagnetic 
interaction certainly present in the cuprates. A possi- 
ble model is the Nearly Antiferromagnetic Fermi Liquid 
model (NAFFL) of Pines and coworkers>ii*i^ It needs 
to be pointed out, however, that there is no a priory 
reason why the electron-phonon interaction itself could 
not lead to a d-wave gap and there exists recent work 
on this possibilityiiii^iiSiiSiiLiSiiS In any case, when d- 
wave symmetry is involved, the spectral function acting 
in the gap channel Eq. IjBlajl need not be the same as 
the one that determines the renormalizations in the uj- 
channel Eq. ljBlb|) . At Tc in the normal state it is only 
the latter spectral density that enters. There exists con- 
siderable literature on extensions of Eliashberg theory to 
include a d-wave gap based on model spectral densities 
for the electron-boson interaction that may be involvedi^ 
Of course, there is no guarantee that the final theory of 
strongly correlated systems that is needed to describe the 
oxides will fall within the class of boson exchange mod- 
els. Nevertheless, such an approach has proven valuable 
in providing insight into the physics of the oxides as we 
will see also in this paper. 

In the recent literature, tunneling spectroscop y^^i^^i^^ 
as well as angular resolved photo emissionSii2S42Si2L2& has 
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been used to analyze data in terms of boson structure. 
Here we wish to concentrate on optical datai22i22i2ii^ 
Optimally doped YBa2Cu3 06.95 (YBCOe.gs) was first 
considered within a complete Eliashberg formalism gen- 
eralized to include c?-wave pairing by Carbotte et al?^ 
(CSB). A model form for the electron-boson spectral 
function coming possibly from exchange of spin fluctu- 
ations and denoted by /^x(u;) is assumed with two fit- 
ting parameters, the coupling and the spin fluctuation 
energy w^j in the model of Millis et alM (MMP) which 
are varied to get the best fit to the normal state infrared 
data at or close to T = Tc. For the superconducting 
state the same form of /^x(^) is assumed to also deter- 
mine the gap channel but its magnitude is different and 
is fit to get the measured value of Tc. In addition, it 
is found that the data in the superconducting state in- 
dicates the formation of an optical resonance in I'^xi^) 
not present at Tc which increases in amplitude as T is 
reduced and is positioned at 41 meV. Similar optical res- 
onances were later found in the superconducting state 
of other cuprates although not in alli^i In some the res- 
onance seems to persist even in the normal statei^ii^^ 
While the work described above involves a least squares 
fit of an assumed form for Pxi^) to the optical scat- 
tering rate data other inversion techniques^ have been 
considered but so far these are based on approximate an- 
alytic formulas for the relationship between the optical 
scattering rate and the electron-boson spectral density 
rather than the full Eliashberg formulation of Carbotte 
and coworkers iiSS*22i2£ 

Such approximate formulas were given by AUen^ for an 
electron-phonon system and are based on ordinary sec- 
ond order perturbation theory at zero temperature. Allen 
considered the normal as well as the superconducting 
state with s-wave symmetry. A generalization to finite 
temperature was provided by Shulga et aL— who only 
considered the normal state but started directly from an 
Eliashberg formalism and the Kubo formula for the con- 
ductivity. A generalization to include as well a pseudo- 
gap was recently provided by Sharapov and Carbotte*^ 
Finally, Carbotte and Schachinger^^ generalized the orig- 
inal work of Allen to a superconductor with d-wave gap 
symmetry. 

The advantage of these simplified but approximate 
equations is that they relate directly through an integral 
the optical scattering rate to the desired spectral function 
I'^xi^) ^-iid various numerical techniques such as singu- 
lar value decomposition's*^ can be used to numerically 
invert the equation. Of course, the alternate method of 
assuming some characteristic functional form for /^x([j) 
and least squares fit a few parameters to the data can 
also be employed based on the simplified equations de- 
scribed above instead of employing the full Eliashberg 
equations. For instance, the new equations of Sharapov 
and Carbotte^ have already been used in this way by 
Hwang et alm^ to analyze data in underdoped Ortho-II 
YBCOe.s- 

The aim of this paper is to understand better how limi- 



tations in accuracy of the simplified formulas can impact 
on the resulting form of I'^xi'^) ^nd to explore as well 
the advantages and limitations of numerical inversion 
techniques such as SVD and Maximum Entropy Method 
(MaxEnt) as well as least squares fit. 

The paper is organized as follows. In Sec.^]the formal 
background is discussed. One subsection concentrates on 
the three major methods of inversion, namely the second 
derivative method, deconvolution methods based on ap- 
proximate relations, and the least squares fit method. 
The second subsection discusses approximate formulas 
for the normal and superconducting state which allow 
one to calculate the optical scattering rate T^p^{uj) from 
a given spectral density a'^F{uj) using a convolution in- 
tegral. Section ITTll discusses numerically the caveats and 
merits of the various methods of inversion by studying 
normal metals as well as High- Tc cuprates. Computer 
generated and experimental t~j^(lu) data for the normal 
and superconducting state are used as input for the inver- 
sion. Finally, conclusions are drawn in Sec. IIVI Two ap- 
pendices have been added. Appendix^gives an overview 
of the Maximum Entropy method in terms of Bayesian 
probability theory. Appendix presents all important 
equations which allow to calculate the optical scattering 
rate within the framework of full Eliashberg theory. 



II. FORMALISM 

A. Methods of Inversion 

In order to understand the mechanism of superconduc- 
tivity it is important to have detailed knowledge of the 
spectral function a'^F{uj) and as tunneling, the estab- 
lished source of information on a'^F{uj) ^ was initially not 
a successful tool in the high- Tc superconductors, the in- 
frared optical conductivity, aop{i-o), became increasingly 
important, particularly in the form of the optical scat- 
tering rate 

r-p\cj) = ^^e[a-p\u)] (1) 

of extended Drude theory. Here, flp is the plasma fre- 
quency. 

There are, in principle, three methods to extract the 
information on a'^F{uj) from the optical scattering rate 
(inversion). An essential requirement for a solution ob- 
tained with any of these methods is that the result should 
match the data points as well as possible. In order to 
assess the quality of the fit we need to know how the 
experimental data points ti = T~J'{uJi) scatter around 
the 'true' values t^ = T~p[ijji), that is, we need to 
know in terms of Bayesian probability theory the like- 
lihood p(t|t°,I). It describes the distribution of N data 
points t = {ti\i — 1,...,7V} given the 'exact' values 
t° = {t,^|i = 1, . . . ,iV}, which are usually expressed in 
terms of the parameters of the physical model. The sym- 
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bol X designates all additionally available background in- 
formation comprising the experimental setup as well as 
the physical model employed. 

The likelihood is determined by the experimental setup 
and we have to keep in mind that the experimental signal 
contains at least three contributions^ 

T-}{uj) = T-p'iuj) + B{lo) ± (2) 

with B(ijj) usually a slowly varying background signal 
which is typical of the experimental setup and r]{Lo), the 
noise in the data. 

Unfortunately, we do not have any knowledge concern- 
ing the functional form of the likelihood for the experi- 
mental data sets considered in this study. Therefore, we 
make the assumption of an uncorrelated normal distri- 
bution with standard deviation a: 

p{t\t\l)^{2na^r^l^e^v{-\l')- (3) 

Here 7^ is the misfit 

l'-T.'l (4) 

i 

which is expressed in terms of the residuals 



They measure the deviation of the data points ti from the 
'true' values t° (or the best estimates thereof) in units of 
the error bar a. 

For lack of further information, we argue that the like- 
lihood lO is a reasonable choice and we note that it is, 
in fact, the most uninformative probability distribution 
given only the mean value and the variance and no fur- 
ther information.'^^ We want to stress that curve fitting 
with minimization of the misfit 7^ is also implicitly based 
on the assumption of a Gaussian likelihood. 

The problem of obtaining a^F{i^) from t~J-{uj) is ex- 
tremely ill conditioned. This implies that a direct solu- 
tion will be totally dominated by noise and, therefore, be 
completely meaningless. For this reason, all methods dis- 
cussed here involve a nuisance or regularization param- 
eter that can be tuned in order to suppress noise contri- 
butions. Apart from ad-hoc settings, a sensible choice is 
to adjust the regularization parameter such that j'^ — N 
is obtained. (See Appendix 1X1) 

The first method of inversion is based on the 
relationshipiS 

which is approximately equal to a^F^uj) in the normal 
state at zero temperature. Application of this formula 
to experimental data will result in numerical difficulties 
because we have to keep in mind that the experimental 



signal T~J-{uj) consists, according to Eq. of at least 
three contributions. Two of these, namely B(uj) and ri(uj) 
can obscure completely the looked for spectral function 
a^F((jj) when the second derivative of t^J-{uj) is calcu- 
lated. 

On first sight, Eq. jnj would require that t~J-{uj) must 
be ambiguously smoothed 'by hand'^^ which is certainly 
not true. First of all, it is much better to 'smooth' 
the function u}t~J^(uj) which is monotonically increasing, 
much less structured, and equal to zero at w = 0. The 
application of standard data processing techniques like 
Fast-Fourier- Transform (FFT) smoothing or EFT low 
pass filters on this function allows one to remove quite 
reliably the noise contribution 77(0;). For instance, the 
upper frequency threshold applied to the FFT low pass 
filter will play the role of a nuisance (or renormaliza- 
tion) parameter in this particular case. If there is further 
knowledge about the background function B(uj) applica- 
tion of Eq. lO is much safer than it looks on first sight. 
We will discuss caveats and merits of this second deriva- 
tive method later on using computer generated results 
which ensure B{u!) = and which allow a controlled 
noise contribution rj^Lo). 

The second method of inversion is based on the decon- 
volution of the approximate relation 

00 

T-/(w;T) = JdnK{uj,n;T)a''F{uj), (7) 


where T denotes the temperature. The kernel K{u!, fl; T) 
is determined from theory. The caveat of this method is 
that solutions of Eq. iQ for a^Fiuj) are not unique be- 
cause, generally, the deconvolution of Eq. {7)) constitutes 
an ill-posed problem. 

There are two approaches to solve this deconvolution 
problem and both are based on a discretization of Eq. {Tj) 
of the form 

N2 

T-j,\u;,;T) = Y,^^jK{uj,,n,;T)a'F{n,), (8) 

with i — 1, . . . , TVi and AQj = Qj-^-i — flj. The first ap- 
proach is straight forward and is called Singular Value 
Decomposition"^^ (SVD) which is based on the vector 
form of Eq. (|SJ| , namely 

t = Ka, (9) 

with the vector t — {ti = T~p^^^{uji; T)|i = 1, . . . , A^i}, the 
matrix K = {Kij = Arjji^'(wi, fi^ ; r)|i = 1, . . .iVi, j = 
1,...7V2}, and the vector a = {aj — a'^F{nj)\j — 
1,...,N2}. Using SVD, the matrix K of dimension 
iVi X N2 is transformed into the matrix product USV , 
with U and V being unitary matrices of dimension 
A^i X and N2 x N2 respectively. The matrix S — 
diag{sj|j = 1, . . . , with Sj the singular values (svs). 
Finally, V'^ denotes the transposed matrix V. If the 
vector t and the matrix K are known the vector a and. 
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thus, a'^F{(jj) can be determined by 'inverting' Eq. 
a = VS'U^twith S' = diag{l/sj|j = 1, . . . , N2}. How- 
ever, noise contained in the data t will be dramatically 
magnified by smallest svs, rendering the result meaning- 
less. For this reason, all contributions by svs below a 
certain threshold have to be discarded by replacing the 
corresponding diagonal elements 1/si in the matrix S' 
by zeros. This threshold plays the role of the nuisance 
parameter in the SVD method. Dordevic et alm^ stud- 
ied this approach extensively and discussed in particu- 
lar the number of svs necessary to get a 'smooth' spec- 
tral function a^F^uj) together with a reasonable recon- 
struction of the input data. In principle the problem of 
'smoothing by hand' is moved from the input to the out- 
put of the process. The caveat of this approach is the fact 
that it doesn't ensure that the resulting spectral function 
a'^F{uj) be positive definite. Most of the time, a^Fiuj) 
will contain negative parts which cannot be removed even 
by applying further regularization schemes^ Such neg- 
ative parts are unphysical. 

The second approach to the deconvolution problem 
is the so-called Maximum Entropy Method (MaxEnt). 
Originally, E.T. Jaynes-^' suggested the Maximum En- 
tropy principle for the assignment of probability distribu- 
tions: If only some testable information such as the mean 
value is given, one should select that probability dis- 
tribution {pi} which maximizes the Shannon entropy^!'' 
S = —^i^iPi^og{pi) subject to all known constraints. 
In the case where only the mean and the variance are 
known, the normal distribution is the 'most uninforma- 
tive' probability distribution (pdf). Although the 'true' 
pdf may be completely different, a normal distribution 
can be a sensible choice for lack of further background 
information. 

The MaxEnt principle has been generalized to the in- 
ference of strictly positive functions such as the spec- 
tral function a'^F{u}) within Bayesian probability theory. 
This fully probabilistic description allows for an explicit 
treatment of the ambiguity inherent in badly conditioned 
problems and is discussed in some detail in Appendix fXl 
In our particular case the generalized Shannon-Jaynes 
entropy (|A3p is applicable with the default model vector 
m chosen to be constant. Most of the time this constant 
is chosen in such a way that the spectral function a^Fluj) 
develops a certain high energy behavior. 

The third method of inversion uses model spectral 
functions which depend on a few parameters which are 
then determined using a least squares fit to experiment 
based either on approximate formulas of the form Q or 
full Eliashberg theory. Very often preliminary results de- 
rived with the help of the second derivative method from 
experiment (or using one of the other above mentioned 
methods) can be utilized to minimize the number of pa- 
rameters to be fitted. Results from other experiments, 
for instance inelastic neutron scattering etc., can easily 
be incorporated. Nevertheless, in general this method 
will also result in non-unique solutions for a'^F{uj). 



B. Approximate Formulas 

For the normal state at zero temperature AUen^ pro- 
vided a simplified form of the kernel of Eq. lO , namely 
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{uj-n)9{uj-n), (10) 



where 9(x) is the step function. This formula is based on 
a second order perturbation theory approach based on 
the weak electron-phonon coupling in normal metals and 
it is valid only in the clean limit, i.e.: no impurities. To 
overcome the zero temperature restriction Shulga et a/.— 
started from a full Eliashberg description of the electron- 
phonon formalism and applied a series of approximations 
to reduce the full results to the approximate form 
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2cjcoth ( — ] - {lu + ri)coth 



2T 



2T 



+ (^ - r!)coth (^^^^) 



(11) 



which properly reduces to Eq. H10() for T = 0. When 
applied to invert data one has to keep in mind that this 
kernel becomes singular for f2 = 0. 

The work of Shulga et al. was generalized recently 
by Sharapov and Carbotte-^ to treat the possibility of 
a pseudogap opening up in the fully dressed density of 
states N{bj). They obtain 



K{u:,n;T) 



TT 



de 



NiO) 



N{0) 



X [n{n;T) + fin-e-T)] [f{e~u;:T) 
-f{e + u;;T)], (12) 

which properly reduces to the result 1)11(1 of Shulga et 
alm^ when N{uj) is taken to be constant and equal to 
N{0). Here n{u!;T) and f{u!;T) are the Bose and Fermi 
distributions, respectively. The zero temperature limit of 
Eq. ((T^ was obtained by Mitrovic and Fiorucci^ based 
on Allen's second order perturbation theory approach. 

Allen also provided a kernel similar to Eq. (|10|1 which 
applies approximately in the superconducting state at 
zero temperature. In this case the kernel is of the form 
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(13) 



It ensures that T^p^{uj) is zero for < w < 2Ao. Here, 
E{x) is the complete elliptic integral of the second kind 
and Ao is the energy gap at T = 0. It is valid in the 
clean limit only. To derive Eq. (|13|l Allen treated the 
superconducting transition within the framework of BCS 
theory, i.e.: Eq. H13|) is only valid for s-wave symmetry of 
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the superconducting order parameter. Moreover, Aq is 
an external parameter to Eq. and its value has to be 
determined by other means. Treating the superconduct- 
ing transition within the framework of Eliashberg theory 
will certainly go beyond the possibilities of Eq. H13|l and 
this will have to be kept in mind when Eq. H13|l is ap- 
plied to invert superconducting state optical data of real 
s-wave superconductors which are well known to be ex- 
ceptionally well described by Eliashberg theoryi^ 

A consensus exists that in the high cuprates the 
superconducting order parameter is of d-wave rather 
than s-wave symmetry. Here we follow Carbotte and 
Schachinger^?. and simulate (in a first approximation) the 
effect of d-wave in the Allen formula H13|l for s-wave by 
simply averaging it over a distribution of gaps having d- 
wave symmetry. The result is that Eq. H13|) needs to be 
averaged over the polar angle of the two dimensional 
Cu02 Brillouin zone. This results in the kernel 
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({lu -n)0{Lu + 2Ao{d) -n) 



xE 



(14) 



with (• • • )ij denoting the "iJ-average which can be limited 
to the interval e [0,7r/4] for symmetry reasons. Fur- 
thermore, Ao('!9) = Ao cos(2i?) reflecting the d-wave sym- 
metry of the superconducting order parameter. Eq. H14(l 
ensures that the optical scattering rate is finite in the 
superconducting state for a; > 0. This is in contrast to 
what is observed in s-wave superconductors. 



III. NUMERICAL RESULTS 
A. Normal Metals 

We will study in quite some detail various inversion 
techniques using, as a first material, lead. The electron- 
phonon spectral density a'^F{uj) was derived from tun- 
neling data by McMillan and Rowellf^ This spectrum, 
which is represented in the top frame of Fig. H by a 
gray solid line, has two distinctive peaks which are sep- 
arated from each other by about 4meV. The Debye en- 
ergy Ljjj = 11.2 meV. Optical data for lead was obtained 
by Joyce and Richards^ and later by Farnworth and 
Timuskn^ The extracted a'^Ftr{u}) was found to be in 
remarkable good agreement with earlier tunneling data 
and with the results of direct band structure calculations 
of a^F{LLj) by Tomlinson and Carbotte.^ In this subsec- 
tion only computer generated optical scattering rate data 
based on the various kernels discussed in Sec. Ill Bl and on 
complete Eliashberg equations (Appendix^ will be used 
to study the various inversion techniques. This provides 
us with data free of a background signal B{u!) and with 
controlled noise 77(0;). 

In a first step zero temperature, normal state T~p{uj) 
data are computer generated using kernel (|10|l . We calcu- 
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FIG. 1: Top frame: Inversion of zero temperature, normal 
state optical scattering rate data Toj^{uj) of lead computer 
generated using the kernel HIO^ . The gray solid line indi- 
cates the a^F{Lij) employed to generate the data. The dotted 
line corresponds to the function W{uj) according to Eq. 
SVD(a;) (dashed line) shows the result of an SVD inversion, 
and ME{uj) (solid line) presents the result of a MaxEnt inver- 
sion. Bottom frame: The crosses correspond to the normal- 
ized uncorrelated Gaussian noise rj{ijj)/G which was added to 
the input data for the MaxEnt inversion and r{ijj) (solid line) 
gives the residual of the MaxEnt data reconstruction. 



late the function W{Ld) (dotted line, upper frame Fig. 
using the second derivative method and the agreement 
with the input spectrum (gray solid line) is almost per- 
fect without the need of 'smoothing by hand'. (Only the 
second peak shows oscillations.) Inversion of the input 
data using the SVD method results in the curve SVD(a;) 
(dashed line, upper frame Fig.P). The svs threshold was 
set to 10"'^, i.e.: 87 svs have been used. A few wiggles 
remain in the valley between the two peaks and we see 
oscillations at energies > 10 meV which also go negative. 
Data beyond the Debye energy are irrelevant. 

For the application of the MaxEnt method we add un- 
correlated Gaussian noise of standard deviation a — 
to ensure a controlled error distribution for the com- 



FIG. 2: Temperature dependence of the normal state optical 
scattering rate T^^(tj) of lead. The solid lines correspond to 
results derived using Eliashberg theory and the dotted lines 
correspond to data generated using kernel lllll . Temperatures 
are 0.3 K, 1 K, 10 K, and 50 K. 

puter generated T~p{uj) data. (This, in principle, bi- 
ases the comparison in favor of the second derivative 
and SVD method.) The curve ME(a;) (soHd line, up- 
per frame Fig. ^ presents the result of the MaxEnt in- 
version in which we used optimized preblur [blur-width 
h = 0.05, see Eq. (|A4|) ] and the default model was set to 
ruj = 0.001. ME(a;) underestimates slightly the second 
peak but otherwise shows perfect agreement with the in- 
put spectrum. We also see an additional feature beyond 
the Debye energy which is irrelevant because it reflects 
the default model. The bottom frame of Fig. ^ demon- 
strates the quality of the data reconstruction achieved by 
the MaxEnt method. The crosses symbolize the normal- 
ized noise rj{uj)/a which was added to the optical scat- 
tering rate and the solid line corresponds to the residual 
As we only added noise to the computer generated 
Tqp^(cl;), riuj) should track the normalized noise ri{uj)/a, 
as it does. 

Zero temperature is not a realistic case and we proceed 
to study normal state, finite temperature results. There 
are two options to computer generate T~p{Lij) data: (a) 
kernel (|11|) is applied, or (b) Eliashberg theory (see Ap- 
pendix is used. The superconducting order parame- 
ter is zero in the normal state and the renormalization 
formula ljBrb|) takes on a closed form. Fig. |21 presents 
our results for the temperature dependence of the opti- 
cal scattering rate in lead for four different temperatures, 
namely 0.3 K, IK, 10 K, and 50 K. The results according 
to Eliashberg theory are presented by solid lines, while 
the dotted hues correspond to the results of kernel ((TT|l . 
There are small but distinct differences between the two 
sets of data. 

Fig.inipresents the spectra W{uj), SVD(u;), and ME(cj) 
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FIG. 3: Inversion of finite temperature, normal state r^^(aj) 
data computer generated using kernel IIH . The solid lines 
correspond to the temperature T — 0.3 K, dashed lines to 1 K, 
dotted lines to 10 K, and dash-dotted lines to 50 K. The gray 
solid line represents the a'^F{uj) spectral function applied to 
calculate the optical scattering rate data. Top frame: Second 
derivative method. Center frame: SVD method. Bottom 
frame: MaxEnt method. 



which result from the application of the first two methods 
of inversion discussed in Sec. Ill Al As input we used the 
finite temperature normal state optical scattering rate for 
lead generated using kernel (fTT)l (dashed lines of Fig.|21). 
The top frame presents as a result of the second derivative 
method, the function W{ll!) as defined in Eq. ©. At the 
lowest temperature, T = 0.3 K the input a'^F{uj) (grey 
solid line) is perfectly reproduced (solid line), while at 
T — IK (dashed line) the high energy peak is already 
underestimated. At T = 10 K (dotted line) the method 
is no longer able to resolve the two peak structure, and at 
T = 50 K (dash-dotted line) the method fails completely. 
Nevertheless, it has to be emphasized that no smoothing 
had to be applied to the input data as no artificial noise 
was added. 
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The center frame of Fig.|3]presents the resuhs SVD(ci;) 
of a singular value decomposition of Eq. using ker- 
nel ini). The svs threshold was set to 10 At T = 
0.3 K and 1 K we obtain reasonable agreement with the 
a^F^uj). For energies > 9meV the inversion shows os- 
cillations and SVD(a;) even becomes negative which is 
unphysical. We also see oscillations at low energies and 
between the two peaks. At T = 10 K the SVD method 
still resolves a hint of a two peak structure in contrast to 
the second derivative method. Finally, at T = 50K the 
method fails. 

The bottom frame of Fig.|31presents the results ME(ci;) 
of the MaxEnt deconvolution of Eq. (jHJ. Uncorrelated 
Gaussian noise of tr = lO"'^ was added to the computer 
generated data. For the inversion optimized preblur (see 
Appendix was applied with b = 0.4 for T = 0.3 K, 
b = 0.46 for IK, b ^ 0.89 for T = lOK, and b = 1.91 
for T = 50 K. The default model was set to nij — 0.01. 
The T — 0.3 K inversion (solid line) gives almost perfect 
agreement with the model a^-F(w) spectral function. At 
T = 1 K the high energy peak is underestimated but 
reproduced at the appropriate energy. At T = 10 K, the 
two peak structure is still well resolved, only the second 
peak is underestimated and shifted to lower energies. The 
result is certainly much better than that of the other two 
methods. Finally, at T = 50 K MaxEnt is no longer able 
to resolve the two peak structure. Nevertheless, it is quite 
interesting to note that the area under the dash-dotted 
curve is 3.95 meV which is very close to the area 4.03 meV 
under the original a^F^uj) spectral function. 

Fig.|21demonstrated that T^p^{uj) data computer gener- 
ated from full normal state Eliashberg theory differ from 
the approximate results of kernel (|ll|l . It can also be as- 
sumed that metals will more likely follow the predictions 
of full Eliashberg theory rather than approximate model 
formulas. It is therefore interesting to investigate how 
the inversion on the basis of the approximate kernel Hll|) 
performs when t~j^(uj) data computer generated within 
full Eliashberg theory are used for inversion. The result 
is presented in Fig. ^ which is organized the same way 
as Fig. 13 The top frame demonstrates the application of 
the second derivative formula which shouldn't have any 
problems because this method is not based on approx- 
imate models formulas. Nevertheless, W{uj) is only in 
reasonable agreement with the original a'^F{u!) at low 
temperatures. The low energy peak is over estimated 
and shifted towards higher energies, the valley between 
the peaks is too low, the second peak is positioned at 
the correct energy but its height is over/underestimated. 
At T = 10 K the two peak structure is no longer re- 
solved and at r = 50K the inversion fails. The central 
frame of Fig. 0| presents the function SVD(a;) as a re- 
sult of an SVD inversion. The svs threshold was set at 
10~^. At low temperatures, the peak positions are at the 
proper energies, nevertheless the low energy peak is over- 
estimated, the valley between the peaks underestimated, 
and the high energy peak is too wide. At T = 10 K a two 
peak structure is resolved but both peaks are placed at 




FIG. 4: The same as Fig. 01 but now the Eliashberg theory 
generated normal state, finite temperature t~j^ (uj) data (solid 
lines in Fig. |5J are used as input. 



the wrong energies. At T = 50 K the method fails alto- 
gether. It is typical for this method to show oscillations 
at energies > 9 meV which result in unphysical negative 
contributions even at energies below lud. 

Before applying the MaxEnt inversion uncorrelated 
Gaussian noise of cr = 0.1 was added to the input data. 
For the temperature 0.3 K and 1 K the preblur parameter 
was optimized to 0.54 at higher temperatures no preblur 
was applied. The default model was set to 0.01. At 
low temperatures MaxEnt overestimates the low energy 
peak and/or makes it broader. The valley between the 
peaks is too low, the second peak is resolved reasonably 
well and is at the proper position but underestimated 
in height. At higher temperatures the two peak struc- 
ture is no longer resolved (T = 10 K, dotted line and 
T = 50K, dash-dotted line). Nevertheless, data recon- 
struction is within error bars and this proves that we face 
in this case a deconvolution problem which is particularly 
ill conditioned. 

All this demonstrates quite clearly that the application 
of methods of inversion based on approximate models to 
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FIG. 5; Top frame: The computer generated superconducting 
state optical scattering rate T~p{ijj) for Pb at T = O.OSTc with 
Tc = 7.2K. The sohd line is based on the full Eqs. (lEJ while 
the dashed line is obtained using the simplified kernel (1131 . 
Bottom frame: The gray solid line symbolizes the a'^F{Lo) of 
Pb. The solid line is the result of the inversion of Top{i-tj) 
data computer generated using kernel 1131 (dashed line in 
the top frame of this figure) using the MaxEnt method while 
the dashed line corresponds to an SVD inversion. The dash- 
dotted and dotted lines present equivalent results but now for 
full Eliashberg data (solid line in the top frame of this figure). 



experimental data (represented here by computer gener- 
ated T~p^{uj) data using full Eliashberg theory) can quite 
easily result in deconvoluted spectra a'^F{u}) which will 
be close but not necessarily equal to the real electron- 
phonon spectrum which governs the interaction. In par- 
ticular, the deviations from the gray solid lines in the 
central and bottom frame of Fig. 0] represent the devia- 
tions from the 'real' Q!^F((jj) which are required by the 
approximate kernels to reproduce the input optical scat- 
tering rate data as well as possible. 

We now move on to a discussion of superconducting 
state data. The top frame of Fig.|Slpresents the results for 
the superconducting state optical scattering rate T~p (uj) 



in Pb at T = 0.05Tc with Tc = 7.2 K. The solid line was 
obtained on evaluation of the full Eqs. HB1|I and HB2() 
taken for s-wave symmetry of the superconducting order 
parameter and a Coulomb pseudopotential /i* — 0.1438. 
The dashed line is for comparison and was obtained from 
kernel (|13|l using the electron-phonon spectral density 
a^Fluj) shown by a gray solid line in the bottom frame 
of this figure. The approximate kernel H13|) is evaluated 
with Aq = 1.39 meV, the gap edge predicted by the full 
Eliashberg calculation. 



The bottom frame of Fig.[Slpresents the result of SVD 
as well as MaxEnt inversions based on the approximate 
kernel (|13|l . The dashed line corresponds to the SVD 
inversion (svs threshold was set at 10^^) of the optical 
scattering rate generated using kernel H13|) (dashed line 
in the top frame of Fig. O and Aq = 1.39 meV. As ex- 
pected, the agreement is almost perfect. The dotted line, 
on the other hand, shows the result of an SVD inversion 
of full Eliashberg data (solid line in the top frame). Here, 
the lower, transverse peak centered around 4meV is 
broader and the area under the peak is larger. The same 
holds for the upper, longitudinal phonon peak but the dif- 
ferences are now less pronounced. Beyond ~ 9 meV the 
dotted line becomes negative giving unphysical results. 
These differences are due to the use of the approximate 
kernel (|13|l in the inversion process. 



The MaxEnt inversion was performed by attaching er- 
ror bars of cr = 10^^ to the data and by adding uncor- 
related Gaussian noise of the same a. Furthermore, no 
preblur was applied and the default model was set to 0.1. 
The solid line is the MaxEnt inversion of the data com- 
puter generated with the help of kernel H13|l (dashed line 
in the top frame of Fig. [SJ. Again we achieve perfect 
agreement. At energies > 10.5 me V the function ME (a;) 
levels off at the value 0.1 demonstrating the influence of 
the chosen default model. The dash-dotted curve, on the 
other hand, is based on full Eliashberg theory generated 
input data (solid line in the top frame of this figure). 
Both peaks are now overestimated in their height and 
width and are shifted towards higher energies. Neverthe- 
less, ME(w) never becomes negative which proves that 
there exists a positive definite solution for the deconvo- 
lution problem of Eq. (jT)) . As experimental data are more 
likely to be close to full Eliashberg theory results, the de- 
convolution of Eq. I7J on the basis of kernel H13|l will re- 
sult in an electron-phonon spectral density a'^F{u!) which 
will not agree in all details with the real spectral density 
despite the fact that the input data will be excellently re- 
produced. The only possible check for the validity of the 
deconvoluted electron-phonon spectral density SVD(w) 
or ME(tt') is using it to calculate the optical scattering 
rate based on full Eliashberg theory and compare with 
the data. Such a comparison will then result in neces- 
sary readjustments of the deconvoluted spectrum. 
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B. High-Tc cuprates 

In contrast to the normal metal lead, the high 
cuprates are not likely to be electron-phonon systems, 
they are known to be highly correlated systems. There is 
a class of models used to describe such systems which we 
will refer to as boson exchange models. They have many 
common elements with the electron-phonon case. In par- 
ticular, there exists a well developed literature on the 
Nearly Antifcrromagnetic Fermi Liquid model (NAFFL) 
introduced by Pines and coUaboratorsi^i^ The exchange 
bosons are antifcrromagnetic spin fluctuations as de- 
scribed by Millis et alM (MMP). Within this model the 
Eliashberg equations are retained in a zeroth order ap- 
proximation neglecting possible vertex corrections which 
go beyond Migdal's theorem. The electron-phonon spec- 
tral density a^F(a;) is replaced by the imaginary part 
of the spin susceptibility multiplied by the square of a 
coupling of the spin fluctuations to the charge carriers. 
In general, this interaction is anisotropic and not pinned 
to the Fermi surface.-^ Nevertheless, as a first approxi- 
mation, one can work with a simple interaction spectral 
function I'^xi^) which replaces the a^F(ijj) of Eliashberg 
theory. Carbotte et al^ found that in optimally doped, 
twinned YBCO single crystals the measured normal state 
optical scattering rate t~J{lo), reported by Basov et al^, 
can be well described by a single MMP form: 



Uj/UJSF 



1 + {uj/uosfY 



(15) 



The two parameters, the square of the coupling constant 
/ and the characteristic spin fluctuation energy los f were 
determined from a least squares fit to the data in the 
energy interval < w < 250 meV. The values P — 0.83 
and ujsF = 20meV have been reported by CSB. 

Based on the results for lead we cannot necessarily ex- 
pect inversions of T^p{Lo) measured around T = 100 K 
or even higher will be feasible. To investigate this, nor- 
mal state T^j^{uj) data are computer generated at various 
temperatures, namely T = 1 K, 10 K, 50 K, and 100 K us- 
ing either kernel Hll|l or full Eliashberg theory with an 
I'^xi^) determined by Eq. (|15|) and with the above values 
for the parameters and losf- The results of the inver- 
sion based on the approximate kernel Hill) are discussed 
in Fig. El with T^p^{u;) data generated using kernel l(TT)) 
and, in Fig. [7| with Tgp{uj) generated by full Eliashberg 
theory. 

The top frame of Fig. presents results for W{llj) from 
the second derivative method. At the two lowest tem- 
peratures, namely T = IK (solid line) and 10 K (dashed 
line) the I'^xi^) spectrum (gray solid line) is almost per- 
fectly reproduced. At higher temperatures, namely at 
T = 50K (dotted line) and 100 K (dash-dotted fine) 
the inverted spectrum develops a less pronounced peak 
which is also shifted towards higher energies. In the tail 
(uj > 100 meV) noise develops in the inverted spectrum. 
Nevertheless, in contrast to lead with its narrow two peak 
structure the simple MMP form can easily be inverted 
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FIG. 6: The same as Fig. Elwith the a^F(a;) replaced by the 
/^X(w) defined by Eq. with 7^ = 0.83 and wsf = 20 meV. 
The solid lines corresponds to T = 1 K, the dashed lines to 
10 K, the dotted lines to 50 K, and the dash-dotted lines to 
100 K. 



from optical scattering rate data even at temperatures 
around 100 K. 

The center frame of Fig. El presents the results SVD(a') 
of a singular value decomposition. The svs threshold was 
set at lO^'^ for T = IK and 10 K and was increased to 
10^2 foj. r = 50K and 100 K. The inverted spectrum 
agrees reasonably well with the original spectrum at low 
energies, w < 75 meV. At higher energies significant os- 
cillations occur. We also note for T — 100 K a typical 
splitting of the peak at 20meV into two peaks. This is 
an indication of a particularly ill conditioned inversion 
problem. This phenomenon can be observed for rather 
narrow and fast rising peaks. 

The bottom frame of Fig.jHlpresents the results ME(w) 
of a MaxEnt deconvolution. The error bars on the data 
were determined by cr = 0.15 and uncorrelated Gaussian 
noise of the same a was added. No preblur was applied. 
The default model was set to 0.05. This ensures that 
the spectrum ME(ci;) agrees at high energies with I'^xi^)- 
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FIG. 7: The same as Fig.inibut now Eliashberg theory gener- 
ated optical scattering rates have been used as input for the 



inversion. 



The agreement with the original I'^xi'^) spectral function 
(gray solid line) is excellent up to temperatures of 50 K. 
At 100 K the peak at 20meV is not as well resolved and 
shows tendency towards a splitted peak similar to the 
SVD inversion at the same temperature. Otherwise, the 
agreement is still rather good. 

The results presented in Fig. [T] are very similar to the 
ones shown in Fig. with the difference that SVD and 
MaxEnt (error bars were determined by cr = 0.2, un- 
correlated Gaussian noise of the same a was added, no 
preblur) now overestimate the peak at 20 meV. The SVD 
result for 100 K (dash-dotted line in the center frame of 
Fig. Cj) develops a slight tendency for a splitted peak at 
20meV. All this is the result of the minor differences 
between the optical scattering rates calculated from ker- 
nel (|ll(l and full Eliashberg theory and the top frame of 
Fig. IHl demonstrates how little these differences are. 

We proceed to a study of the MaxEnt inversion of ex- 
perimental data and make use of the T = 95 K normal 
state optical scattering rate measured by Basov et al."^^ 
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FIG. 8: Top frame: Experimental normal state optical scat- 
tering rate {u) at T — 95 K for an optimally doped, 
twinned YBCOe.gs single crystal as reported by Basov et aim— 
(solid gray diamonds). The black solid line corresponds to the 
MaxEnt reconstruction based on kernel IIH using the spec- 
tral function ME(ti^) shown in the middle frame of this figure 
(black solid line). The dashed line presents t~j^{uj) of a full 
Eliashberg calculation based on the spectral function I'^xi^) 
shown by a gray solid line in the middle frame of this figure 
and the dotted line corresponds to t~J'{uj) computer gener- 
ated using kernel Ulljl and the same I x{^)- Bottom frame: 
The residual r(a;) according to Eq. I14I I of the Ma:xEnt data 
reconstruction. 



on an optimally doped, twinned YBCOg.gs single crystal. 
Fig. IHl presents the result of a MaxEnt inversion. We as- 
sume the experimental data t~J-{uj) to be contaminated 
by a substantial uncorrelated Gaussian noise of cr — 2.5. 
The inversion is performed on the basis of kernel Hll|l 
and the resulting spectral function ME(ti;) is shown as 
a solid line in the middle frame of Fig. |S| ME(a;) then 
replaces a^F(cj) in Eq. iQ which is used to calculate the 
reconstructed optical scattering rate shown by a solid 
line in the top frame of Fig. |S1 It reproduces excellently 
the experimental data (gray solid diamonds). For com- 
parison, the top frame of this figure contains two more 
results, namely T~p^{uj) calculated from full Eliashberg 
theory (dashed hue) using the Pxi^) reported by CSB, 
namely P = 0.83 and losf = 20meV (gray solid line in 
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the center frame of Fig.lS)). The dotted hne, on the other 
hand, corresponds to T~p{w) calculated from Eq. iQ us- 
ing kernel Ullfl and the same Obviously, the two 
results are very close with the dotted line slightly above 
the dashed one at low energies. The opposite holds for 
high energies. This is in agreement with the result found 
for lead (see Fig. . 

Finally, the bottom frame of Fig. |S1 shows the residual 
r([j) which is a measure for the quality of the data recon- 
struction. Apart from the low energy region the recon- 
struction is within the assumed standard deviation (in- 
dicated by the two straight lines at 1 and -1) of cr = 2.5. 
The result is to be compared with the r{Lj) shown in the 
bottom frame of Fig.^which results from the reconstruc- 
tion of computer generated data with additional uncor- 
related Gaussian noise. The r{uj) in the bottom frame 
of Fig.|Slis a rather smooth function which contains for 
energies > 70meV very little stochastic elements which 
could be identified as noise. The various data points ap- 
pear to be rather correlated an effect which could either 
be attributed to an additional background function B{uj) 
[see Eq. (jJJ] or to a 'real' signal. Nevertheless, what is 
important here is the fact that the inverted ME(a;) has 
a nonzero contribution even at energies > 150 meV thus 
establishing a high energy background in /^x(a;) as pre- 
dicted by CSB. Finally, as both spectral functions pre- 
sented in the central frame of Fig. |S1 reconstruct the ex- 
perimental data (solid gray diamonds in the top frame of 
Fig. IHl equally well, they can be used as valid spectral 
functions because of the non-uniqueness of the decon- 
volution problem. Further calculations and comparison 
with other experiments than optical conductivity may 
then help to discriminate between these two spectra. It 
is interesting to point out that the area under /^x(cl;) 
(42 meV) is approximately reproduced by the area under 
the spectrum ME((jj) (41.4 meV) which could be used to 
explain the oscillations in ME(w) as a result of the en- 
hanced main peak at 20 meV. 

Tu et alm^ measured the optical scattering rate of op- 
timally doped Bi2Sr2CaCu208+5 (Bi2212) single crystals 
at various temperatures. They derived, using data anal- 
ysis different from the methods discussed here, that even 
in the normal state at 100 K a resonance peak is seen in 
the function W{oj) while it is rather featureless at 295 K. 
Schachinger and Carbotte^ also analyzed these data us- 
ing a combination of the second derivative method and 
least squares fits to the data. In particular, they found 
that the T = 295 K data are well described by an MMP 
form H15|l in the energy region < w < 250 meV. The 
least squares fit determined the parameters = 0.655 
and UJSF — 82meV using full Eliashberg theory in the 
fitting procedure. 

As MaxEnt turned out to be a rather powerful inver- 
sion technique we revisit the Bi2212 data analysis. We 
assume the experimental data of Tu et al. to contain 
uncorrelated Gaussian noise of cr = 2.0 for T — 100 K 
and cr = 3.0 for T = 200 K and 295 K. The preblur pa- 
rameter was set to 5 for T = 100 K and to 10 for the 
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FIG. 9: Top frame: Experimental optical scattering rate 
data^^ T~^(uj) of Bi2212 (gray solid symbols) reconstructed 
using the MaxEnt method (black solid line) for various tem- 
peratures, namely T = 100 K, 200 K, and 295 K. The dotted 
lines correspond to data generated by full Eliashberg the- 
ory using the J^x('^) spectra reported by Schachinger and 
Carbotte.''^ Bottom frame: The spectral function ME(aj) 
(black solid line for 295 K, black dashed, black dash-dotted 
for 200 K, and black dotted for 100 K) as a result of the Max- 
Ent inversion of the experimental data (gray solid symbols in 
the top frame of this figure). The gray lines (solid for 295 K, 
dotted for 200 K, and dashed for 100 K) show the /^x('^) spec- 
tra reported by Schachinger and Carbotte.^^ The black dash- 
dotted hne presents results of a MaxEnt inversion of the 200 K 
data (solid gray triangles in the top frame) using the I^xiio) 
(gray dotted line) as the default model. 



other two temperatures. The default model was set to 
0.1. The inversion is based on the application of the ap- 
proximate kernel (|ll|l . Fig. |51 discusses the results of our 
calculations. The top frame presents the data reconstruc- 
tion and the bottom frame the inverted spectral function 
ME(a;) in comparison to /^x(a;) spectral functions sug- 
gested by Schachinger and CarbottenS^ It is quite clear 
that the data reconstruction (black solid lines in the top 
frame of Fig.|51) is in excellent agreement with the original 
data (gray solid symbols) at all temperatures. The black 
dotted lines correspond to [lS) data generated from 
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full Eliashberg theory using the spectral func- 

tions presented in the bottom frame of Fig. |51 by gray 
lines, namely solid for 100 K, dashed for 200 K, and dot- 
ted for 295 K. The full Eliashberg results follow the data 
rather nicely in the energy range < uj < 250 meV and 
then deviate systematically to smaller values for energies 
> 250 meV. Thus, the MMP form alone is not sufficient 
to explain the energy dependence of t~J-{lu) in the whole 
energy range < w < 400 meV. 

Comparing the spectral functions ME(ijj) to the /^x(ijj) 
spectra demonstrates rather good agreement at low en- 
ergies for T = 100 K (black dotted line, gray solid line). 
Both spectra show a pronounced peak at 43 meV and 
they even agree in height and width of the peak which 
is rather fortuitous. At T — 200 K (black dashed line, 
gray dashed line) both spectra develop a less pronounced 
peak with the peak in ME(a;) shifted away from 43meV 
to higher energies. Such a shift towards higher energies 
with increasing temperatures has already been observed 
in the analysis of computer generated r^^(w) data for 
YBC06.95, Fig. El and this peak in the ME{uj) could very 
well correspond to a 43 meV peak in the 'real' spectrum. 
Finally, at T = 295 K (black solid line, gray dotted line) 
both spectra agree in showing a rather flat MMP like 
structure peaked around 82 meV with no indication of a 
resonance peak. This analysis corroborates the results re- 
ported by Tu et alS^ and by Schachinger and Carbottei^ 

It is quite important to notice that all ME(a;) spectra 
develop a second structure of comparable height around 
~ 260 meV for all temperatures. Such a structure is miss- 
ing in the Pxi^) spectra. This additional structure is 
required for a faithful reconstruction of the experimen- 
tal T~J-(uj) data at higher energies. To include one more 
check we repeated the inversion of the T = 200 K r^^^ (lu) 
data using MaxEnt but now without preblur and with 
the default model set to the I'^xi'^) spectrum for 200 K 
(gray dashed line in the bottom frame of Fig. ^ . The 
result is presented by the black dash-dotted line in the 
bottom frame of Fig. O In this case the optical reso- 
nance can be found at 43 meV in contrast to the cal- 
culation with the constant default model (black dashed 
line) but it is now significantly enhanced. Even the kink 
in the I'^xi^) spectrum at about 60meV is reproduced 
in ME{uj). What is important, though, is the fact that 
the additional high energy structure around 260 meV ap- 
pears again with approximately the same strength as in 
all other results. Thus, it seems to be a real and new fea- 
ture not captured by a simple MMP form and this proves 
that the charge carrier-exchange boson spectral function 
I'^xi^) ill the cuprates has non-zero contributions up to 
at least 400 meV, a property which cannot be explained 
by a pure phonon mechanism. 

It is interesting to note in closing that, for instance, 
ME(ijj) for T = 295 K can be described using a simple 
model, namely the original MMP form which is replaced 
for u! > 170 meV by a second MMP form peaked at 
U!SF = 260 meV. A least squares fit to t~^(w,T = 295 K) 
using full Eliashberg theory provides an — 0.55 for 
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FIG. 10: The residual r{uj) of the data reconstruction pre- 
sented in the top frame of Fig. |5| The top frame is for 
T = 100 K, the middle frame for 200 K, and the bottom frame 
for 295 K. 



this second MMP form and agreement with the data is 
achieved over the whole energy range < w < 400 meV. 
This is a second, independent proof of the existence of 
this additional high energy structure in I'^x{ll>) for opti- 
mally doped Bi2212. 

In Fig. ^] the residual r{u!) of our analysis of the 
Bi2212 data is presented. The top frame is for 100 K, the 
middle frame for 200 K, and the bottom frame for 295 K. 
The residual clearly shows a stochastic component which 
is much smaller than the assumed values for a in the en- 
ergy region 100 < w < 350 meV. It can be identified as a 
noise contribution. There is obviously, another slowly os- 
cillating contribution to r{uj) which is almost identical in 
its frequency dependence at 100 and 200 K but it doubles 
its period at 295 K. This contribution is very likely to be 
a background signal B(uj) generated by the experimental 
equipment. 

We proceed to investigate the inversion of supercon- 
ducting state T~p{ui) data for YBCOe.gs. We use for 
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computer generated r^^ (w) data the spectral function 
/^x(w) reported by CSB for T = 10 K. It was derived 
from experimental superconducting state r^^(w) data re- 
ported by Basov et alm^ a.t T — 10 K for an optimally 
doped, twinned YBCOe.gs single crystal. This Px{^) is 
based on the normal state Pxi^) for YBCOg.gs which is 
an MMP form ^ with P = 0.83 and ujsf = 20meV 
(gray squares in Figs. El and d). Superimposed is a pro- 
nounced peak at = 41 meV which was found by ap- 
plying the second order derivative method to the exper- 
imental data. The final form, shown using gray solid 
squares in the middle and bottom frame of Fig. ^2 was 
established by a fit of full Eliashberg rj^^ (ui) results to ex- 
periment. This spectrum which extends to 400 meV was 
applied by Schachinger et alM' within full Eliashberg for- 
malism. The authors demonstrated that it was possible 
to reproduce numerous superconducting state properties 
of optimally doped YBCOg.gs in their temperature and 
energy dependence within experimental errors. The need 
for /^x(tj) to extend to several hundred mcV has also 
been reviewed by Basov and Timusk.^^ 

The top frame of Fig. ^] presents computer generated 
superconducting state t~p^{llj) data at T = 10 K as a 
function of energy. The solid line gives the result of a full 
Eliashberg calculation using the solutions of Eqs. ljBl|) 
and ljB2|l on the basis of the spectral function I^xi'^) 
just described. Eliashberg theory also provides a value 
for the zero temperature gap amplitude Aq = 22.03 meV. 
The dashed line presents the optical scattering rate as 
calculated using kernel 1)14(1 . the above value for Aq, and 
the same spectral function Pxi^)- The two results differ 
substantially in the energy region 70 < w < 200 meV. 

The results of an SVD inversion are shown in the mid- 
dle frame of Fig. ^2 The solid line presents the spectral 
function SVD(cj) found from inverting the full Eliashberg 
results (solid line in the top frame of this figure) on the 
basis of kernel (|14|l using Aq = 22.03 meV as an external 
parameter. The dashed line corresponds to the inversion 
of the scattering rate generated by kernel (|14|l using the 
same value for Aq (dashed line in the top frame of this 
figure). In both cases the svs threshold was set to 10^^. 
The agreement of both spectra SVD(a;) with the original 
/^x(a;) (gray solid line) is rather poor keeping in mind 
that the inversion is based on computer generated data. 

The bottom frame of Fig.^Jis organized as the middle 
frame of this figure. It presents spectra ME(a;) as a result 
of a MaxEnt inversion. For the inversion of t~p^{lu) rep- 
resented by the dashed line in the top frame of this figure 
an error bar of a = 0.01 was assumed and no noise was 
added to the data. The inversion was performed using 
historical MaxEnt with 7^ ~ Ni as convergence criterion. 
The resulting spectral function ME{lu) is represented by 
a dashed line. The agreement with the spectrum I'^xi^) 
(solid gray fine) employed to generate T~p^{uj) is perfect as 
was to be expected. The inversion of full Eliashberg the- 
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ory generated r^, 



data (solid line in the top frame of 



this figure) based on the same /^x(w) is less successful as 
the resulting spectrum ME((jj) (solid line) demonstrates. 
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FIG. 11: Top frame: The computer generated superconduct- 
ing state optical scattering rate r^^(a;) at T = 10 K calcu- 
lated from full Eliashberg theory (solid line) or from kernel 
(1141 (dashed line) using the I'^xi'^) spectrum shown as the 
gray solid lines in the middle and bottom frame of this fig- 
ure. Middle frame: The gray solid line corresponds to the 
spectral function I^xi^) reported by CSB for superconduct- 
ing YBCO6.95 at r = 10 K. The solid line shows the spectral 
function SVD(tj) as a result of an SVD inversion of the full 
Eliashberg result (solid line in the top frame of this figure). 
The dashed line shows the same but now the scattering rate 
generated by kernel (1141 is used as input. Bottom frame: The 
same as the middle frame. ME{u)) is the result of a MaxEnt 
inversion. 



For the inversion an error of cr = 0.7 was attached to the 
T~p {uj) data but we did not add noise. Furthermore, Aq 
had to be reduced to 21 meV in order to keep the peak 
in ME(a;) at 41 meV. The peak height is now grossly un- 
derestimated and the spectrum ME{uj) does no longer 
reproduce the normal state background spectrum. Such 
a result was expected because of the pronounced differ- 
ences particularly in this energy region between the full 
Eliashberg theory generated data and the data generated 
using kernel (|14ll . 

We proceed and study the application of the MaxEnt 
inversion based on kernel H14|) using experimental su- 
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FIG. 12: Top frame: The superconducting state optical scat- 
tering rate t~J'{uj) at T = 10 K for an optimally doped, 
twinned YBCOe.gs single crystal. The solid gray diamonds 
present the original data by Basov et aZi—. The solid line 
shows the MaxEnt reconstruction of the input data while 
the dashed line shows the result of a full Eliashberg theory 
calculation based on the spectrum Px{^) reported by CSB 
for T = 10 K. Bottom frame: The gray solid line presents 
the I^x{^) spectrum suggested by CSB for superconducting 
YBC06.95 at T = 10 K. The solid line shows the spectrum 
ME(aj) as a result of the inversion of the experimental data 
by Basov et al}-' represented by solid gray diamonds in the 
top frame of this figure. 



perconducting state t^^{uj) data. The top frame of of 
Fig. El presents the original data by Basov et ai^- re- 
ported for an optimally doped, twinned YBCOg.gs single 
crystal at T = 10 K (solid gray diamonds). The inver- 
sion of this data is performed using historical MaxEnt 
based on kernel (|14|l with 7^ = Ni as criterion of con- 
vergence. An error bar of cr = 3.5 was attached to the 
data and the default model was set to 0.05. The result- 
ing spectrum ME(cj) is presented as black solid line in 
the bottom frame of Fig. 1121 This spectrum was found 
using Ao = 21 meV. This allowed to place the main peak 
in ME(ijj) at 41 meV. It is obvious that the inverted spec- 
trum ME(tj) differs substantially from the /^x(a;) spec- 
trum (gray solid line) reported by CSB. It shows peak 
splitting of the resonance peak at 41 meV and the low 
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FIG. 13: The same as Fig. 1121 but now for a superconducting 
optimally doped Bi2212 single crystal at T = 6K. The data 
have been reported by Tu et al.^-^ The additional dash-dotted 
line in the top frame describes the result of a full Eliashberg 
calculation based on the ME(tj) spectrum which was down 
scaled by a factor 0.95. 



energy peak (w < 5meV) in ME(a;) is caused by an at- 
tempt to use MaxEnt to extrapolate to very small ener- 
gies which are not supported by the r^"^ (uj) data. Nev- 
ertheless, the data reconstruction (black solid line in the 
upper frame of Fig. I12II is excellent. We added for com- 
parison (dashed line in the upper frame of Fig. [T^ the op- 
tical scattering rate as generated by full Eliashberg the- 
ory using /^x(ci;). The agreement with the data does not 
seem to be good enough to justify the particular shape of 
Pxi^) discussed above. Nevertheless, one has to keep in 
mind that all calculations presented here are performed 
without including impurities, i.e.: in the pure case limit. 
Adding impurities improves the agreement between full 
Eliashberg theory and experiment substantially^ 

As a last example we present the reconstruction of ex- 
perimental superconducting state t~J-{uj) data reported 
by Tu et al^ for an optimally doped Bi2212 single crys- 
tal at T = 6K. The results are presented in Fig. [TSI 
which is organized the same way as Fig. 1121 For the 
MaxEnt data reconstruction an error bar determined by 



15 



a = 3.0 was attached to the data. Historical MaxEnt 
with 7^ = Ni as criterion for convergence was appUed. 
The default model was set to 0.2. The inversion is based 
on kernel 1)14(1 and the spectrum I'^xi'^) reported by 
Schachinger and Carbotte'^'^ for T = 6K is shown as a 
gray solid line in the bottom frame of Fig. ll3l for compari- 
son. It contains a peak at 43 meV and an MMP form \ib\ 
as background with P = 0.655 and ujsf — 82meV. In 
this case the agreement between the inverted spectrum 
ME(w) (solid line in the bottom frame of Fig. I13|l and 
Px{^) is much better in comparison to YBCOe.95. This 
confirms the analysis of Schachinger and Carbotte^S, as 
well as a previous analysis of Bi2212 data by Schachinger 
and Carbotte^ based on data published by Puchkov 
et al.'^ We added one more result to the top frame of 
Fig. ^1 presented by a dash-dotted line. It corresponds 
to the result of a full Eliashberg calculation based on the 
ME(ijj) spectrum (down scaled by a factor of 0.95) in- 
stead of /^x(w). The agreement with experiment is now 
very good. In particular, the 'overshoot' right after the 
main rise in the optical scattering rate^ is better resolved 
than in the original Eliashberg calculation (dashed line). 
Thus, the existence of this overshoot is an indication that 
a dip exists in the /^x(<^) spectrum immediately follow- 
ing the resonance peak. This dip will, of course, not 
be as pronounced as it appears in ME(u;) because the 
major part of it stems from the differences between the 
approximate kernel (|14|l and full Eliashberg theory. If 
the ME(a;) spectrum is, furthermore, employed to calcu- 
late the zero temperature gap amplitude Aq within full 
Eliashberg theory a value of 32 meV is found, in excellent 
agreement with experimental resultsiSLS We also note 
that the optical scattering rate t~^^ (w) in the top frame of 
Fig. 1121 (gray solid diamonds) develops a moderate over- 
shoot following the main rise around 80 < w < 110 meV. 
Thus, also in this case not all of the dip which follows the 
resonance peak in the ME(tj) spectrum can be attributed 
to the differences in the approximate kernel ((14|l and full 
Eliashberg theory. 



C. The least squares fit method 

It has been pointed out in the previous subsection 
that the least squares fit method has already been ap- 
plied rather successfully to invert /^x(a;) spectra from 
experiment using full Eliashberg theory together with 
additional information gathered by other means. This 
method is rather clumsy to handle and time consuming 
as one cannot develop a closed algorithm which allows 
one to fit parameters directly given some standard de- 
viation (T which plays the role of the nuisance parame- 
ter. Therefore, we want to study the least squares fit 
method based on the approximate kernel (|14|l for the 
superconducting state of a d-wave superconductor. The 
T~p^ (w) data are generated from full Eliashberg theory for 
the superconducting state at T = lOK using for Pxi^) 
an MMP form lO with P = 0.83 and ujsf = 20meV 
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FIG. 14: Top frame: The residual r{uj). The solid line corre- 
sponds to a least square fit of an MMP form to computer gen- 
erated optical scattering rate data generated by full Eliash- 
berg theory, Aq — 24.03 meV and a = 2.0 was assumed. 
The dashed line corresponds to a MaxEnt analysis of the 
same data for a — 1.75 keeping Aq fixed and the dotted 
line is for a = 0.05 and Aq — 15 meV. The main part of 
the various r{uj) curves is within ±1 indicating perfect data 
reconstruction within the assumed error. Bottom frame: The 
spectrum I'^xi'^) is^sty solid line) is the input spectrum for 
the full Eliashberg calculation. The solid line gives the in- 
verted spectrum as a result of a least squares fit to the data, 
the dashed line the result of a MaxEnt inversion both with 
Ao = 24.03 meV. Finally, the dotted line shows the spectrum 
ME(aj) which resulted from an optimal data reproduction us- 
ing MaxEnt but now Aq = 15meV. 



(gray solid line in the bottom frame of Fig. 114(1 . The zero 
temperature gap Aq = 24.03 meV. The least squares fit 
method is now applied to determine P and ujsf of an 
MMP form by a least squares fit to T~p^{uj) in the en- 
ergy region < oj < 250 meV. The error bar attached 
to the input data is given by ct = 2.0. (This particular 
value of the standard deviation appears to be a realistic 
value for the reconstruction of experimental data as was 
demonstrated in the previous subsection.) No noise was 
added. A consistent data reconstruction was achieved 
with the parameters P = 1.03 and ujsf = 14meV. This 
becomes apparent from Fig. E| in which the results of 
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the least squares method are illustrated. The solid line 
in the top frame of this figure shows the residual r{uj), 
Eq. which is on average well within the assumed a. 
This insures a correct data reconstruction. The bottom 
frame of this figure compares the least squares fit spec- 
trum (solid line) to the original I'^xi^) indicated by a 
gray solid line. It is interesting to compare the areas un- 
der these two spectra, they are 41.5 meV and 42meV, 
respectively, a difference of about 1%. The parameter A 
which is two times the first inverse moment of I'^xi^) is 
also a good parameter to compare. We get A — 2.48 and 
3.09, respectively. 

Figure. ^] contains additional information. We use 
the MaxEnt method to generate an 'educated guess' for a 
later least squares fit to data based on full Eliashberg the- 
ory. If we use Aq = 24.03 meV and assume a = 1.75 his- 
torical MaxEnt reproduces the input data almost equally 
well as our least squares fit. (Dashed line in the top frame 
of Fig. El) The resulting spectrum ME(ijj) (dashed line 
in the bottom frame of Fig. ll4() has its peak at a slightly 
higher energy (~ 24meV) as compared to the original 
/^x(a;) but otherwise, the input spectrum is reproduced 
rather well, albeit not by an MMP form. The area un- 
der this spectrum is 40meV and A = 3.22, again close 
to the result of the least squares fit 'inversion'. We also 
include, for comparison, the result of a MaxEnt deconvo- 
lution with the emphasis on optimal data reconstruction. 
We reduce the error bar on the input data to a = 0.05 
and use Aq as a parameter to be adjusted in order to 
achieve this goal. An almost perfect reproduction is pos- 
sible if Ao is reduced to 15meV. This becomes apparent 
from the residual r{uj) shown as a dotted line in the top 
frame of Fig.^l The resulting spectrum ME(a;) is pre- 
sented by a dotted line in the bottom frame. The peak 
is now shifted to much higher energies, it is wider, and is 
greater in height as compared to the original I'^xi^)- The 
area under the spectrum is 40.2 meV and A = 1.92. This 
demonstrates clearly the approximate nature of kernel 
(|14|l and that a value for Aq found from optimum data 
reproduction is not physically significant. Instead, one 
has to treat Aq as an external parameter to the inversion 
process which is to be determined by other means like, 
for instance, scanning tunneling microscopy . ^ ^ i 

While the various methods of inversion described in 
this section lead to significant differences in the value of 
A obtained, the area under the spectral density varies 
very little. The differences in A can be traced mainly 
to variations in the position of the main fluctuation peak 
and, therefore, it is recommended that independent infor- 
mation on the position of this peak be used, for example 
from neutron scattering. 

As a result of this study one can say that given addi- 
tional information, like the value of the zero temperature 
gap amplitude, both methods, least squares fit and Max- 
Ent, result in comparable spectra, nevertheless, they dif- 
fer qualitatively and quantitatively from the 'real' spec- 
trum /^x(u;). This emphasizes the role of additional in- 
formation beyond the optical data for a successful data 



analysis. 



IV. CONCLUSION 

There exists a well established formalism that relates 
the electron-phonon spectral density a'^F{uj) to the in- 
frared conductivity. It applies to the superconducting as 
well as normal state and involves the Eliashberg equa- 
tions plus a Kubo formula which gives the optical con- 
ductivity CTop(w) from Green's functions. While such a 
formalism is not as well justified in the case of other 
boson exchange mechanisms such as spin fiuctuations, 
it has, nevertheless, been useful to apply it as a first 
approximation with appropriate essential modifications 
such as d-wave gap symmetry. The resulting equations 
are, however, rather complicated and simplified, lowest 
order perturbation theory expressions for the relationship 
between spectral density J^x(w) and optical conductivity 
have played an important role particularly if a main aim 
is to extract qualitative rather than quantitative infor- 
mation on the size and main features of the /■^x(<^) for a 
given set of optical data. If, however, accurate quantita- 
tive information is desired a full Eliashberg formulation 
cannot be avoided. In this paper we provided compar- 
ison between numerical results for the optical scatter- 
ing rate t^j^{lu) based on the exact equations and results 
generated from several often used approximate relations 
between conductivity and spectral density including a re- 
cent generalization which applies to a superconductor at 
T = with d-wave symmetry. 

Another important issue discussed in detail is the ac- 
curacy, advantages, and limitations of various numerical 
methods which are needed to invert data even within the 
limitations of approximate formulas for the optical scat- 
tering rate. These equations relate the optical scatter- 
ing rate measured in infrared experiments to the desired 
spectral density /^x(tj) through an convolution integral 
involving a known, specified kernel K{fl, lu; T) multiplied 
by l'^x{^)- A second derivative technique of the optical 
scattering rate, often used, is also considered. 

In the normal state at low temperatures the second 
derivative technique applied to computer generated op- 
tical scattering rates based on Eliashberg theory repro- 
duces well the spectral function a^F(tj) of Pb. There is 
a slight overestimate of the longitudinal and transverse 
phonon peaks as well as a small shift to higher energy. 
The region between the peaks is slightly underestimated 
and unphysical tails occur beyond the Debye energy but 
these are to be ignored. If the same data is inverted us- 
ing either SVD or MaxEnt, the agreement between the 
spectral function of Pb and and the deconvoluted spectral 
functions remains good even though both inversion meth- 
ods are based on an approximate lowest order pertur- 
bation theory expression for the relation between spec- 
tral function and scattering rate while the T~p{uj) data is 
based on Eliashberg relations. It is noted, however, that 
MaxEnt does somewhat better than SVD which intro- 
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duces additional oscillations into the resulting spectral 
density not present in the a'^F{u}) of Pb. As the temper- 
ature is increased in the normal state all three methods 
begin to fail and at T = 50 K the two peak structure of 
the Pb a^F(ijj) can no longer be resolved. 

For the cuprates in their normal state an often used 
spectral function /^x(cj) is the MMP form of the NAFFL 
model. This function is characterized by a spin fluctu- 
ation frequency losf and a coupling constant between 
spin susceptibility and the charge carriers. The MMP 
form has a peak at a; = losf but is rather smooth and 
extends to energies of order 100 meV. For such a rela- 
tively unstructured spectrum, even at 100 K, the decon- 
voluted spectral function is much closer to the original 
MMP form than was the case for Pb. For the SVD in- 
version, however, there is a spurious splitting of the peak 
at ui — ujsF (for small ujsf) and some additional oscilla- 
tions occur. These oscillations are also seen in the case 
of the MaxEnt inversion but they are less prominent and 
quite small. Of the two inversion methods considered 
above MaxEnt is to be preferred because it has the ad- 
vantage that few assumptions, namely a default model 
and error bars attached to the data, will result in a de- 
convoluted spectrum which allows excellent data recon- 
struction. This spectrum can also be expected to contain, 
at least qualitatively, all the main features of the 'real' 
spectrum. 

As a result of all this, inversion of normal state t^J-{llj) 
data, even at high temperatures, can lead to useful quan- 
titative results if the boson exchange spectral function is 
expected to be rather smooth with little structure. The 
application of methods based on the approximate kernel 
(|ll|l is to be favored because there are only little differ- 
ences between this kernel and full Eliashberg theory. 

For the low temperature superconducting state use of 
the Eliashberg equations for the calculation of the op- 
tical scattering rate T~p{u!) along with mmicrical inver- 
sion based on approximate lowest order perturbation the- 
ory simplified formulas leads to larger quantitative differ- 
ences between the a'^F{uj) spectrum of lead and decon- 
voluted spectrum than were noted for the normal state 
of lead. Nevertheless, both SVD and MaxEnt methods 
yield very useful, qualitative and even semi quantitative 
information on the shape and size of a'^F{oj) with the 
MaxEnt method, again, to be preferred. 

Inversion of experimental data on the optical scatter- 
ing rate (as opposed to computer generated data) already 
exists in the literature in a few cases and these are based 
on full Eliashberg solutions and the Kubo formula. These 
inversions, however, proceed through a least squares fit 
of the scattering rates assuming a specific mathemati- 
cal form for I'^xi^) characterized by a few parameters, 
namely an MMP form with a low frequency cutoff plus 
a resonance peak at a specified frequency. Such fits have 
had considerable success when applied to normal and su- 
perconducting state in optimally doped cuprates. At- 
tempts have been made since additional complications 
arise due to the emergence of the pseudogap. As yet no 



consensus exists as to the origin of this pseudogap and, 
thus, modeling its effect remains controversial. 

Even though, as noted above, the SVD and MaxEnt 
methods are limited due to inaccuracies in the approxi- 
mate second order perturbation form used to relate opti- 
cal scattering rate to spectral density (instead of the com- 
plete Eliashberg analysis), we have used MaxEnt inver- 
sions to confirm the previously obtained least squares fits. 
All qualitative features of the resulting Px{^) function, 
namely, coupling to an optical resonance at low energy 
and to a background extending even beyond 400 meV are 
confirmed and new features have been unveiled. This 
demonstrates the usefulness of such techniques which 
make no a priory assumption on the shape or size of the 
underlying spectral density. To get best results the value 
of the zero temperature gap amplitude Aq (obtained in 
some other way) should be used rather than varied arbi- 
trarily to get a best fit. The numerical spectra for /^x(ijj) 
could be used as a first step in a more complicated in- 
version process which would involve further refinements 
along the lines of the least squares fit procedure described 
above. Other constraints such as known properties of 
the superconducting state could be added as well in the 
fit. It is not clear, however, that this is necessary and 
worthwhile for the oxides where the mechanism is not 
the electron-phonon interaction and additional compli- 
cations such as a breakdown of the Migdal theorem may 
arise. 
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APPENDIX A: THE MAXIMUM ENTROPY 
METHOD FOR DATA ANALYSIS 

The direct inversion of Eq. ((SJ constitutes an ill-posed 
problem. Therefore, there are many different 'solutions' 
(varying orders of magnitude) that fit the data within 
the error bars. The most general solution to this prob- 
lem is the calculation of the posterior probability distri- 
bution (pdf) p{a\t,T) of possible solutions a given the 
data t and all additionally available background infor- 
mation X, i.e.: the matrix K which is defined by the 
underlying theoretical model, the background function 
B{lo), the noise contribution ?/(aj), etc. Bayesian proba- 
bility theory**^ provides the consistent framework for such 
a fully probabilistic description. Bayes' theorem provides 
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the relation 



p(a|t,X) = 



p(t|a,T)p(a|T) 
p{t\J) ■ 



(Al) 



which relates the posterior p(a|t,X) to the likelihood pdf 
p(t|a,X) and the prior pdf p(a|X). Finally, the denomina- 
tor p(t|J) ensures proper normahzation of the posterior. 
The likelihood comprises the model definition and the 
error statistics of the data. Its knowledge is an essen- 
tial prerequisite for any data analysis. The prior, on the 
other hand, should incorporate all available information 
of the problem at hand. In the particular case discussed 
here, the only known constraint is the positivity of the 
function values aj = a^F{Uj). 

Skilling^*^ showed that the most uninformative prior in 
this case is the Maximum Entropy prior: 



N2 



p{a.\a,2) — exp(Q;S') j JJ^ 



(A2) 



Here, a is a renormalization (nuisance) parameter and S 
is the generalized Shannon-Jaynes entropyi^ 
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It measures the distance of the candidate vector a 
from the so-called default model vector m — {mj\j — 
1, . . . , A^2}, which represents the most probable solution 
prior the observation of any data. In case of insufficient 
background information it should be chosen constant, 
i.e.: rrij = const, Vj. Nevertheless, it is adamant to check 
its influence on the solution, as certain features of the so- 
lution might not be supported by the data but instead 
just reflect the initial assumption of the default model. 

The regularization parameter a determines the rela- 
tive influence of the prior compared to the likelihood. In 
the limit a — > 00 one obtains a ^ m as the most prob- 
able solution; for a — + 0, on the other hand, one gets 
the maximum likelihood solution which will be mean- 
ingless for ill-conditioned problems. Within conventional 
approaches, regularization parameters such as a are of- 
ten fixed by hand. Apart from ad-hoc settings, a sensible 
choice is to adjust the regularization parameter such that 
the expectation value of the misfit 7^ is reproducedi^ In 
case of an A^-dimensional uncorrelated normal distribu- 
tion the misfit is described by the ^^-distribution with 
N degrees of freedom and has mean (j^) = N and vari- 
ance var(7^) = 2N. Historically, the criterion -f^ — N 
was employed first in order to fix the parameter a (his- 
torical MaxEnt). However, one has to keep in mind that 
the solution might change dramatically if the regulariza- 
tion parameter a is tuned such that 7^ varies between 
N - V2N < 7^ < TV + V2N. 

In principle, the regularization parameter a can be de- 
termined consistently within Bayesian probability theory 



by computing the most probable value a which maxi- 
mizes the probability ip(a|t, T) g iven the data t. (This is 
the classical MaxEnt of R.ef. 152.) Unfortunately, the cal- 
culation of p(a|t,X) involves high dimensional integrals 
which can only be evaluated using rather crude simpli- 
fications. The approximation usually applied^ tends to 
overfit the data as p(a|t,I) is systematically overesti- 
mated for small a which results in a too small d- value at 
which p{a\t,X) has its maximum as a function of a. Von 
der LindenS^ suggested a different approximation scheme 
that partly corrects these deficiencies and yields results 
similar to the historic criterion. 

For some data sets analyzed in the study, we found 
that all methods to determine the value of a suffered 
from oscillations ('ringing') due to overfitting. This has 
been observed for other applications as welL.^*'^^ To a cer- 
tain extent, this ringing is intrinsic to the MaxEnt prior 
which explicitly treats all points of the reconstruction a 
as uncorrelated. 

In order to enforce smoothness of the solution 
Skilling^^ suggested the introduction of a 'hidden image' 
h which is blurred by a Gaussian 
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Here, the Xj designate the abscissas of aj and hj. The 
vector a enters the likelihood, while h is used to compute 
the entropy 5*. The blur-width b is an additional hyper- 
parameter that can be determined simultaneously with 
a by locating the maximum of p{b, a\t,X) given the data 
t in the spirit of Ref. [5^ 

Various choices of the blur-width b can be regarded 
as distinct models which have a different number of de- 
grees of freedom (similar to fit functions involving differ- 
ent numbers of parameters). For 6 ^ all positive dis- 
crete representations a can be realized as a — s- h, while 
in the limit b 00 only constant functions Oi = const, 
can be represented, i.e.: the model has only one effective 
degree of freedom. 

The optimal blur-width b is determined by the inter- 
play of the likelihood and Occam's razoriM*^ which pe- 
nalizes the complexity of the model employed and is im- 
plicit in the calculation of p{b, a\t,2). The 'penalty fac- 
tor' is the ratio of the width of the likelihood and the 
prior distributions. Thus, a simpler model may be more 
favorable because a larger fraction of the parameter space 
is likely to be realized according to data although a more 
complex model fits the data better. 

Unless stated otherwise, we have determined the op- 
timal blur-width b for the MaxEnt reconstructions pre- 
sented in Sec. IIIII as outlined above. For the computa- 
tion of p(6|t,X) we chose a flat prior p(6|I) on the interval 

bmin <b< bmax with bmin ~ X2-X1 and bmax ~ Xn-Xi- 

The MaxEnt method obviously allows for an explicit 
treatment of ambiguous solutions and it allows prior 
knowledge to be taken into account consistently by intro- 
ducing a suitable prior pdf. A direct inversion, like the 
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SVD method, which may be badly conditioned or may APPENDIX B: ELIASHBERG EQUATIONS 

involve uncontrolled approximations, is avoided. Finally, 
it is possible to obtain error estimates. Nevertheless, it 
has to be pointed out that 'fuzzy' constraints such as 
smoothness of the output of the inversion process make 
a definition of the prior pdf rather complicated^ 

I 



The generalized, clean limit Eliashberg Equations which play an important role in this study are 

oo 

A{v + iO+;e) ^ nTg ^ cos{29) [X{iy - iujjn) + + ii^,n)]h{iuJm) (Bla) 



m=0 
oo 

+iTTg / dz cos(2^)/^x(z) [n{z) + f{z — i^)] h{iujrn ^ v — z -\- iO'*'), 



and, in the renormalization channel, 

oo 



m=0 

oo 

+in J dz l'^x{z)[n{z) + f{z - v)]g{iuj„, ^ - z + iO+). (Bib) 

— oo 



Here 



, A{iujm; 0) cos(2e) \ ,■ \ w(ia;,„) 



and the parameter g allows for a possible difference in spectral density between uj and A channels. It is fixed to 
get the measured value of the critical temperature. In the above /S.(iu}m]0) is the pairing energy evaluated at the 
fermionic Matsubara frequencies LVm — TrT(2m — l),m = 0, ±1, ±2,...; f{z) and n(z) are the Fermi and Bose 
distribution, respectively. The renormalized Matsubara frequencies are ULi{iLUm)- The analytic continuation to real 
frequencies ly of the above is A(i^ + iO~^; 9) and uj{oj + iO~^). The brackets {■ ■ ■)s are the angular average over 9, and 
X{iy) — J^^dV,a^F{n)/{i' ~V, + iO^). Eqs. (|B1|I are a set of nonlinear coupled equations for the renormalized pairing 

potential A(i/ + iO+; 9) and the normalized frequencies ui{i^ + iO+) with the gap A(i/ + iO^; 9) — /S.{v + iO+; 9)/Z{v), 
where the renormalization function Ziy) was introduced in the usual way as Cj{v + zO+) = vZ(y). To get the s-wave 
version of these equations g is set equal to one and all cos(20) factors are to be omitted with no average over the polar 
angle 9. A Coulomb pseudopotential /i* must also be introduced in Eq. IjBlap . 

The optical conductivity follows from knowledge of Co and A. The formula to be evaluated is 

aop(r, JdLo tanh ( l^) [J{u;, v) - J{-lo, ) . (B2) 



2 



The function J(w, i') is given by 

2J{uj,v) 



1 - N{lu; 9)N(lo + u;e)- Pjco; 9)P{lu + iv; 9) 
E{u-9)+E{uo + v-9) 
1 + N*{uj; 9)N{uj + u]9)+ P*{uj; 9)P{uj + 
^ E*{uj]9) - E{u + v;9) 



(B3) 



with E{uj]9) = Y^w2(w + i0+) - A2(tj + iO+;6'), N{u]9) = Cj(bj + iQ+)/E{bj;9), and P{ijJ]9) = h{uj + iQ+;9)lE{uj;9). 
Finally, the star refers to the complex conjugate. 
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